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Abstract 

From the first stages of star and galaxy formation, non-gravitational processes such 
as ram pressure stripping, SNs, galactic winds, AGNs, galaxy-galaxy mergers, etc... 
lead to the enrichment of the IGM in stars, metals as well as dust, via the ejection 
of galactic material into the IGM. We know now that these processes shape, side 
by side with gravitation, the formation and the evolution of structures. We present 
here hydrodynamic simulations of structure formation implementing the effect of 
the cooling by dust on large scale structure formation. We focus on the scale of 
galaxy clusters and study the statistical properties of clusters. Here we present our 
results on the Tx — M and the Lx — M scaling relations which exhibit changes 
on both the slope and normalization when adding cooling by dust to the standard 
radiative cooling model. For example, the normalization of the Tx — M relation 
changes only by a maximum of 2% at M = 10^^ whereas the normalization of 
the Lx — Tx changes by as much as 10% at Tx = 1 keV for models that including 
dust cooling. Our study shows that the dust is an added non-gravitational process 
that contributes shaping the thermodynamical state of the hot ICM gas. 
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1 Introduction 

The high m etallicity observed in t he intra-cluster medium (ICM hereafter - see 
for instance iBorgani et al.l (120081 ). for recent works) is understood within ana- 
lytical models and numerical simulations as a consequence of the various grav- 
itational (galaxy-galaxy mergers, galaxy interactions, ram pressure stripping) 
and non-gravitational (SN powered galactic winds, AGN outbursts, starburts 
winds) processes at play in the ICM. These processes drive the formation and 
the evolution of galaxies within their environment (see for insta nce works by 



Kapferer et all (120061 ) : lDomainko et al.l (120061 ) : lMoll et al.l (120071 ). Moreover, it 
is obvious that the process of tearing material from galaxies does not discrim- 
inate between metals, gas, stars and dust. All are injected in the ICM/IGM 
further on leading to the metal enrichment of the surrounding medium, and 
possibly to its dust enrichment. To date the searches for a diffuse dust compo- 
nent within the ICM have not been significantly successful. Toget her with their 



study of the reddening of background galaxies towards clusters, iMuUer et al 



( 120081 ) gathered the existing constraints on the presence of dust in the ICM 
and the derived constrain ts in terms of extinction or dust-to-gas mass ratio. 
Montier and Giard (l2005h obtained a clear statistical detection for the total 
dust emission in clusters (from member galaxies and ICM) by using the stack- 
ing of the IRAS data (at 60 and 100 /im near the peak of dust emission) at 
the lo cation of more tha n 11000 clusters and groups. In an extension of this 
work, I Giard et al. (l2008h have investigated for the first time the evolution of 
the total IR bolometric luminosity of clusters with redshift. Comparing to 
the expected IR luminosity from the integrated SFR in cluster galaxies and 
considering the low end of current c onstraints for the d ust-to-gas mass ratio 
in the ICM (IChelouche et al.l . l2007h . ICiard et all (l2008h concluded that ICM 



dust could represent 2% of the bolometric IR luminosity of clusters. On the 
other hand, from a theoretical point of view the ef fect of dust on a ICM/IGM - 
type thermalized plasma has been formalised by iMontier and GiardI (120041 1. 
These authors have computed the cooling function of dust within an optically 
thin hot thermal plasma taking into account the energetic budget for dust. 
In view of their results, dust thus comes within the ICM/IGM, as an added 
non-gravitational process that may influence the formation and evolution of 
large scale structures in a non-negligible way. 

In order to address this problem, we have put into place the first hydrodynamic 
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Fig. 1. Total cooling (radiative + dust) functions implemented in the numerical 
simulations (black lines). The radiative and dust cooling functions alone are re- 
spectively shown in blue and red. The figure illustrate the dependency of the dust 
cooling efficiency as a function of the metallicity if the IGM (red curves on both 
panels) and as a function of the dust-to-gas mass ratio within the ICM (left and 
right panels). 

A^-body numerical simulations of hierarchical structure formation implement- 
ing the effect of dust cooling according to the dust nature and abundance. In 
this proceedings paper, we present original preliminary results of this study 
focussing on the properties of clusters of galaxies at redshift zero. 



2 Cooling the ICM with dust 



Montier and GiardI (120041 ) have shown that in the temperature range of galaxy 
clusters, dust can be a cooling agent of the ICM, which can compete efficiently 
with radiative cooling. The strength of cooling by dust depends on dust grain 
properties and more specifically on dust abundance and dust grain sizes. The 
more abundant the dust and the smaller the grains, the higher the efficiency 
of the cooling by dust is. 

• The dust abundance, i.e the dust-to-gas mass ratio, is a function of the rel- 
ative abundance of dust with respect to the solar vicinity dust abundance 
(i.e Zci^q) and of the metallicity of the ICM (i.e Z/Zq) : Zd = Mdust/Mgas = 
fd (Z/Zq) Zd^Q- We investigate three values for fd- 0.001, 0.01 and 0.1. The 
corresponding dust-to-gas mass ratios bracket roughly the current theoret- 



10.-^ and 10-3 ( 


Popescu et al.. 


200d; 


Aguirre et al.. 


2001; 


Chelouche et al.. 


2007; 


Muller et al.. 


2008|; 


Giard et al. 


20081). 



At a given metallicity, the most relevant parameter for the cooling function 
by dust is the grain size. Indeed, the smaller the grains, the higher the 
cooling power of the dust is. We chose to test three types of size distribution: 
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Fig. 2. Density- Temperature phase space diagram of the radiative cooUng simulation 
(left) and a case of radiative+dust cooling simulation. The separation of the different 
phases (ICM, IGM, cool gas) is illustrated by the dashed-yellow lines. The overlaid 
contours are isoprobability contours of 0.025, 0.05, 0.1, 0.15 and 0.2. 

two fixed grain sizes with a = 10~^ /im and a = 0.5 /im, respectively labelled 
small and big. The third assume s for th e IGM dust grains a distribution in 

-3.5 



sizes as defined by iMathis et al.l (119771 ) for the galactic dust: N{a) cx a 



within the size interval of [0.001, 0.5] /xm. It is further referred as the 'MNR' 
distribution. 

Thus in our simulations, once the metallicity is known, a and fd are the only 
two parameters driving the dust cooling rate (i.e A(a, Z) = A{a,fd)). Fig. 1 
illustrates the dependency of the dust cooling function with respect to the 
metallicity and dust abundances for a MNR dust grain size distribution. 



3 Simulation of structure formation with dust cooling 



3.1 Hydrodynamic N-body simulations 

We have run a set of seven numerical simulations using the public code 
package Hydra, an adaptive particle-particle /particle -mesh (AP^M), gravity 



solver with smoothed particle hydrodynamics (S PH) fjCouchman et al.l . Il995 



Couchman . 1991 : Pearce and Couchmanl . 1997). The global gas metallicity 
is assumed to evolve with time as Z = 0.3(t/to)ZQ, where Zq is the solar 
metallicity and t/to is the age of the universe in units of the present time. All 
simulations are issued from the same initial condition snapshot, at z = 49. The 
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initial density field was constructed, using = 4, 096, 000 particles of bary- 
onic and dark matter, perturbed from a regular grid of fixed comoving size 
L = 100 h~^Mpc. We assumed a A-CDM cosmology with parameters, Q = 0.3, 
Qa = 0.7, Qb = 0.0486, h = 0.7. The amplitude of the matter power spec- 
trum was normalised using as = 0.9. With this choice of parameters, the dark 
matter and baryon particle masses are 2.1 x 10^° /i~^Mq and 2.6 x 10^ h~^MQ 
respectively. The gravitational softening in physical coordinates was 25 h~\pc 
below z = 1 and above this redshift scaled as 50(1 + z)~^ h~^kpc. Gas in the 
simulations is allowed to cool using the total table presented in Fig 1. At a 
given time step, gas particles with overdensities (relative to the critical den- 
sity) larger than 10^, and temperatures below 1.2 x lO^K are converted into 
collisionless baryonic matter and no longer participate in the gas dynamical 
processes. 

3.2 Phase space and cluster maps 

On Fig. 2 we present the phase space diagram for the whole box at 2; = for 
the radiative cooling (left panel) and the radiative-dust cooling model with 
a = 0.001/im and fd = 0.1 (right panel - model Dl in Table 1). 

The displayed quantity is the probability distribution function of gas parti- 
cles (computed as the number of particles in logarithmic bins of size 0.1 dex 
divided by the total number of particles) as a function of temperature and 
overdensity. The colour scale is linear, as represented on the top of the figure. 
The features in the diagrams reflect the gas physical effects included in the 
simulations: gravitational heating, adiabatic expansion, radiative cooling and 
cooling by dust (the latter effect is only present on the right panel). Note that 
our simulations do not include energy feedback nor photoinization-heating due 
to a UV background. The absence of the later effect causes the low density, 
low temperature, gas to cool adiabatically without being regulated by pho- 
toionization heating. For T > 10^ K gas cools more efficiently according to 
the cooling rates of Fig. 1. At high overdensities log(p/ < p >) ^ 2.5 gas is 
mainly distributed in bound objects such as galaxy groups and clusters with 
high temperatures. The gas in these regions cools rapidly due to the absence 
of energy feedback. The inclined V-shape (in black) on the far right end side 
of the diagrams is a region of fast cooling, where the gas cannot remain over 
a Hubble time. In both panels, the two upper peaks at intermediate overden- 
sities reflect the shape of the cooling functions in Fig. 1. The contribution 
of cooling due to metals {Z/Zq = 0.3 at 2; = in our simulations) increases 
the cooling efficiency at T = 10^ K as compared to a gas with primordial 
composition. The absence of non-gravitational heating in our simulations also 
enhances the appearance of these two upper peaks in the phase space diagram. 
The horizontal peak at T ~ 10^ K arises because in our simulations the gas 
below this temperature can only cool adiabatically. At high overdensities, gas 
that cools below 10^ K is considered collisionless baryonic matter and it is 
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Fig. 3. Projected density (top pane) and temperature (bottom) maps for one of the 
most massive clusters found in the simulations. The radiative case is presented on 
the left hand side, whilst a model implementing radiative+dust cooling is presented 
on the right hand side of the figure. Density and temperature maps are respectively 
in logarithmic and linear scales. The overlaid circle's radius is the virial radius of 
the cluster. The X and Y axis are given in arcmin considering a projection of the 
clusrter positionned at z = 0.03 (in this case the spatial scales are just indicative, 
but meaningless). 

converted into galaxy fragments/ star forming material. 

The differences between the two diagrams are not striking. However, the con- 
centration of gas in the cool and low-density phase of the Dl model is higher 
than in the purely gas cooling case (model C). For the ICM gas, the phase 
space distributions are also different when looked in detail. More particles are 
present in the form of cold gas (i.e bellow the horizontal yellow lines) in the 
radiative+dust cooling direct consequence of the extra cooling from 

the dust component. As a consequence, the hot and less dense gas distribution 
(typically located in the clusters core, i.e. top-right parts of the diagrams) is 
more spread towards slightly less hot and dense gas in the radiative+dust 
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Table 1 

Parameters of the simulated models presented on Fig 4. 



Run 


r liybluS 


Id 


vjriaiii size 


c 


cooling (no dust) 






Dl 


cooling with dust 


0.100 


small 


D2 


cooling with dust 


0.100 


MRN 


D4 


cooling with dust 


0.010 


MRN 



cooling case. This is due to a well-known effect related with the condensation 
of gas that is converted into coUisionless star forming material in the simula- 
tions. The removal of cold gas from the hot phase simply leads to an increase 
of the mean local temperature of the remaining gas. 

For the same models we display in Fig. 3 the projected density and tempera- 
ture maps for one of the most massive cluster found in our simulations. As for 
the phase space, the changes in features on the density maps are not obvious. 
However, the locus of high density seems brighter in the radiative+dust cool- 
ing case, e.g the cluster centre is brighter and more peaked. The temperature 
maps present more visible differences between the dust /no-dust case. Some 
hot areas are enhanced or have appeared when adding dust cooling. This is 
also a consequence of the aforementioned (counter intuitive) effect of addi- 
tional cooling (e.g due to dust) actually contributes to raise the temperature 
of the remaining ICM gas. 



4 Cluster scaling relations 

4-1 Tx — M and Lx — M at redshift zero 

In fact, to be able to properly quantify the effect of the cooling by dust in the 
ICM, we have investigated the statistical properties of the cluster populations 
forming in our simulations. In each simulation we selected halos with a total 
mass larger than 5 x 10^^ M0. We ended up with a minimum of 60 clusters 
in each simulations at 2; = 0. Here we illustrate our study with two scaling 
relations given at z = 0: the (emission weighted) temperature-mass relation, 
Tx — M, and the X-ray bolometric luminosity-mass relation, Lx — M. The 
two relations are presented on Fig. 4 for the radiative cooling case, and three 
radiative+dust models described in Table 1. 

The Tx — M relation is the less disturbed of the two. Results for model D4 
is almost compatible with the radiative case. Model D2 and Dl illustrate the 
effects of the dust grain size distribution and the dust abundance. 
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Fig. 4. (left) Emission weighted temperature-mass relation, Tx — M, (right) X-ray 
bolometric luminosity-mass relation, Lx — M, at redshift zero. Symbols and best 
fit lines are color-coded for each presented model (see Table 1). The grey shaded 
area corresponds to the ilo" dispersion for the radiative model (i.e model C). 



These small variations are in fact understandable because the temperature of 
clusters is essentially dominated by gravity and because the total cluster mass 
(that is dominated by the dark mater component) is in practice unaffected by 
the extra condensation of a fraction of the cluster gas mass. However small, 
we note that the emission-weighted vs mass scali ng is more sens i tive t o dust 
than the mass- weighted vs mass relation studied in Ida Silva et al.l (120091 ). This 
confirms expectations, because Tx gives more weighting to gas with stronger 
X-ray emission, i.e gas with high density and temperature locate in the clus- 
ter central regions. As a consequence the Tx — M scaling not only exhibits 
larger changes of amplitude but also larger differences in the slopes because of 
the relative impact of dust on Tx for different mass systems. Altogether the 
variation in normalisation are less than ^ 2% over the 3 presented models at 
M = 10^^ Mq. 

On the other hand the Lx — M relation is typically a baryon driven scaling 
law. Which means that it is very sensitive to the physics of baryons at play 
in the ICM. Again, going from model D4 to Dl via D2, one can witness the 
effect of the grain size (as for Dl and D2 the abundances are the same), and 
then the abundance effect from model D2 to D4. Both the abundance and 
the grain size distribution have similar amplitude, and the overall changes in 
normalisation at 1 keV goes up by ~ 10%. It is also striking how the dust 
cooling efficiency increases with the ICM temperature. Taking into account 
the range of investigated temperatures, this is perfectly consistent with the 
dust cooling function imp lemented in the simulations, thus with the study of 



Montier and GiardI (120041 ) that shows this discriminating effect with the tem- 



perature. 

Also, we have pointed out that these changes induced by dust cooling on the 
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normalisation as well as the slope of the scaling relations, are strongly linked 
to the physical properties of dust. We have illustrated that the effect of the 
dust grain size is as efficient as the effect of the abundance of dust grains 
populating the ICM. 



4 .2 Limitations of the model 



It is obvious that our present implementation of dust cooling in simulations 
has intrinsic limitations and should be considered as a first ("zero order") 



approach to the problem, see Ida Silva et all (120091 ) . Indeed, with our implen- 
tation of dust in our simulations, we directly linked the presence (thus abun- 
dance) of dust with the presnce of metals (i.e we put dust wherever there 
are metals in the simulations). However, dust and metals are not sensitive 
to the same physical proc esses, and if the metals are not destroyed, the dust 
grains are sputtered. From lDraine and Salpeterl (119791 ) and lMontier and Giard 
( 12004| ). for a hot thermal plasma, i.e 10^ < T < 10^ K, and densities met in 
the core of dense clusters, i.e uh ~ 10~^ cm~^, and grain sizes ranging from 
0.001 and 0.5 /im, dust grains have lifetimes of 10^ < tdust < 10^ yr. The over- 
all efficiency of cooling by dust in the ICM/IGM is thus strongly linked to the 
physical processes of dust injection. Therefore, beside the cooling function of 
dust, with which we have performed a fully self-consistent implementation of 
the effect of dust as a cooling vector of the ICM/IGM, we did not performed 
stricto senso a physical implementation, as we do not deal with the physics of 
the dust creation and dust destruction processes. This would be the next step 
in our study. 

Also, we directly correlated the abundance of dust with metallicity. The evolu- 
tion law we chose for the metallicity does underestimate the metallicity at high 
redshift. Indeed, Z = 0.3(t/to)ZQ normalised at 0.3^0 at 2 = gives ~ 0.2 at 
z = 0.5 and ~ 0.1 at 2 = 1. However, observations as well as theoretical studies 
on the metal enrichment of the ICM have i ndicate high values of the metal 



icity, i.e Z ~ 0.3Zq, up to redshift above 1 (jCora et al.l . l2008l : iBorgani et al. 



20081 ). An early enrichment of the ICM/IGM with dust, then sputtered, could 
(partly) explain these high values of the metallicity. Therefore, by underes- 
timating the metallicity at high redshifts, we may have underestimated the 
amount of dust injected in the ICM at high redshifts, and thus the efficiency 
of dust cooling when integrated from an early epoch down to redshift zero 

They are probably other limitations to our results, such as the local thermal 
equilib rium hypothesis underlied in the computation of iMontier and Giard 



( 20041) cooling function (following a non-equilibrium hypothesis ( iDwek et al 
1990l )). that may lead to difference in the cooling curves, the used theoretical 



setup for the dust emission within an optically thin hot thermal plasma is well 
suited to our purpose of assessing the influence of dust on the galaxy clusters 
properties. 
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5 Conclusion 



With our current "zero" order implementation of the effect of dust cooling, 
we have demonstrated the actual effect of the presence of dust within the 
IGM/ICM as a non-gravitationnal process that plays a role on the physical 
properties of baryons within large scale structures. Deeper studies including 
other scaling laws such as the entropy-mass, S — M, and the int egrated Comp- 



tonisa tion parameter-mass, Y — M, relations are presented in Ida Silva et al 
fl2009h . 
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